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Abstract 

The new ekpyrotic scenario attempts to solve the singularity problem by involving vio- 
lation of the null energy condition in a model which combines the ekpyrotic/cyclic scenario 
with the ghost condensate theory and the curvaton mechanism of production of adiabatic 
perturbations of metric. The Lagrangian of this theory, as well as of the ghost condensate 
model, contains a term with higher derivatives, which was added to the theory to stabilize 
its vacuum state. We found that this term may affect the dynamics of the cosmological 
evolution. Moreover, after a proper quantization, this term results in the existence of a new 
ghost field with negative energy, which leads to a catastrophic vacuum instability. We explain 
why one cannot treat this dangerous term as a correction valid only at small energies and 
momenta below some UV cut-off, and demonstrate the problems arising when one attempts 
to construct a UV completion of this theory. 
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1 Introduction: Inflation versus Ekpyrosis 

After more than 25 years of its development, inflationary theory gradually becomes a standard 
cosmological paradigm. It solves many difficult cosmological problems and makes several predic- 
tions, which are in a very good agreement with observational data. There were many attempts to 
propose an alternative to inflation. In general, this could be a very healthy tendency. If one of 
these attempts will succeed, it will be of great importance. If none of them are successful, it will 
be an additional demonstration of the advantages of inflationary cosmology. However, since the 
stakes are high, we are witnessing a growing number of premature announcements of success in 
developing an alternative cosmological theory. 

An instructive example is given by the ekpyrotic scenario pQ. The authors of this scenario 
claimed that it can solve all cosmological problems without using the stage of inflation. However, 
the original ekpyrotic scenario did not work. It is sufficient to say that the large mass and entropy 
of the universe remained unexplained, instead of solving the homogeneity problem this scenario 
only made it worse, and instead of the big bang expected in pQ, there was a big crunch [2], [3]. 

Soon after that, the ekpyrotic scenario was replaced by the cyclic scenario, which used an 
infinite number of periods of expansion and contraction of the universe Unfortunately, the 
origin of the scalar field potential required in this model, as well as in [Tj, remains unclear, and the 
very existence of the cycles postulated in [1] has not been demonstrated. When this scenario was 
analyzed using the particular potential given in [I] and taking into account the effect of particle 
production in the early universe, a very different cosmological regime was found [3 [6]. 

The most difficult of the problems facing this scenario is the problem of the cosmological 
singularity. Originally there was a hope that the cosmological singularity problem will be solved 
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in the context of string theory, but despite the attempts of the best experts in string theory, this 
problem remains unsolved [Tl, El E] - Recently there were some developments in the analysis of this 
problem using the AdS/CFT correspondence [10], but the results rely on certain conjectures and 
apply only to five dimensional space. As the authors admit, "precise calculations are currently 
beyond reach" for the physically interesting four dimensional space-time. This issue was previously 
studied in [TTJ , where it was concluded that "In our study of the field theory evolution, we find 
no evidence for a bounce from a big crunch to a big bang." 

In this paper we will discuss the recent development of this theory, called 'the new ekpyrotic 
scenario' [121 EES EE2 US] , which created a new wave of interest in the ekpyrotic/cyclic ideas. This 
is a rather complicated scenario, which attempts to solve the singularity problem by involving 
violation of the null energy condition (NEC) in a model which combines the ekpyrotic scenario 
P] with the ghost condensate theory [16] and the curvaton mechanism of production of adiabatic 
perturbations of metric [T71 IT8]. 

Usually the NEC violation leads to a vacuum instability, but the authors of [121 [131 EEH E5] 
argued that the instability occurs only near the bounce, so it does not have enough time to fully 
develop. The instability is supposed to be dampened by higher derivative terms of the type — (D0) 2 
(the sign is important, see below), which were added to the action of the ghost condensate in |16j . 
These terms are absolutely essential in the new ekpyrotic theory for stabilization of the vacuum 
against the gradient and Jeans instabilities near the bounce. 

However, these terms are quite problematic. Soon after introducing them, the authors of the 
ghost condensate theory, as well as several others, took a step back and argued that these terms 
cannot appear in any consistent theory, that the ghost condensate theory is ultraviolet-incomplete, 
that theories of this type lead to violation of the second law of thermodynamics, allow construction 
of a perpetuum mobile of the 2nd kind, and therefore they are incompatible with basic gravitational 
principles [HI [201 12H 122] ■ 

These arguments did not discourage the authors of the new ekpyrotic theory and those who 
followed it, so we decided to analyze the situation in a more detailed way. First of all, we found 
that the higher derivative terms were only partially taken into account in the investigation of 
perturbations, and were ignored in the investigation of the cosmological evolution in [T2 | \13 \ [HI [15] . 
Therefore the existence of consistent and stable bouncing solutions postulated in the new ekpyrotic 
scenario required an additional investigation. We report the results of this investigation in Section 

El 

More importantly, we found that these additional terms lead to the existence of new ghosts, 
which have not been discussed in the ghost condensate theory and in the new ekpyrotic scenario 
[T2"l [T51 dH [T51 [12]. In order to distinguish these ghosts from the relatively harmless condensed 
ghosts of the ghost condensate theory, we will call them ekpyrotic ghosts, even though, as we will 
show, they are already present in the ghost condensate theory. These ghosts lead to a catastrophic 
vacuum instability, quite independently of the cosmological evolution. In other words, the new 
ekpyrotic scenario, as well as the ghost condensate theory, appears to be physically inconsistent. 
But since the new ekpyrotic scenario, as different from the ghost condensate model, claims to 
solve the fundamental singularity problem by justifying the bounce solution, the existence of the 
ekpyrotic ghosts presents a much more serious problem for the new ekpyrotic scenario with such 
an ambitious goal. We describe this problem in Sections El 0, SI and [71 



2 



Finally, in Appendix we discuss certain attempts to save the new ekpyrotic scenario. One of 
such attempts is to say that this scenario is just an effective field theory which is valid only for 
sufficiently small values of frequencies and momenta. But then, of course, one cannot claim that 
this theory solves the singularity problem until its consistent UV completion with a stable vacuum 
is constructed. For example, we will show that if one simply ignores the higher derivative terms 
for frequencies and momenta above a certain cutoff, then the new ekpyrotic scenario fails to work 
because of the vacuum instability which is even much stronger than the ghost-related instability. 
We will also describe a possible procedure which may provide a consistent UV completion of the 
theory with higher derivative terms of the type +(D0) 2 . Then we explain why this procedure fails 
for the ghost condensate and the new ekpyrotic theory where the sign of the higher derivative 
term must be negative. 



2 Ghost condensate and new ekpyrosis: The basic sce- 
nario 

The full description of the new ekpyrotic scenario is pretty involved. It includes two fields, one of 
which, 0, is responsible for the ekpyrotic collapse, and another one, \i is responsible for generation 
of isocurvature perturbations, which eventually should be converted to adiabatic perturbations. 
Both fields must have quite complicated potentials, which can be found e.g. in [T5]. For the 
purposes of our discussion it is sufficient to consider a simplified model containing only one field, 
0. The simplest version of this scenario can be written as follows: 



1 /D0 X 1 



where X = is dimensionless. P{X) is a dimensionless function which has a minimum at 
I/O. The first two terms in this theory represent the theory of a ghost condensate, the last one 
is the ekpyrotic potential. This potential is very small and very flat at large 0, so for large this 
theory is reduced to the ghost condensate model of [16]. 

The ghost condensate state corresponds to the minimum of P{X). Without loss of generality 
one may assume that this minimum occurs at X = 1/2, i.e. at <9j0 = 0, = — m 2 , so that 
= —m 2 t. As a simplest example, one can consider a function which looks as follows in the 
vicinity of its minimum: 

P(X) = i(X -1/2) 2 . (2) 

The term — \{^r) 2 was added to the Lagrangian in [16] for stabilization of the fluctuations of the 
field in the vicinity of the background solution 0(t) = — m 2 t; more about it later. 

This theory was represented in several different ways in [TBI [T2l [TBI [Til [15], where a set of 
parameters such as K and M = M 2 /M' was introduced. The parameter K can always be absorbed 
in a redefinition of M; in our notation, K = 1. 
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The equation for the homogeneous background can be represented as follows: 



dt 



,x 



+ 



d t (4> + 3Hj>) \ 



a 3 V,0 m 4 
~M~ 4 



where we introduced the notation 



m 



M'M 2 

9 9 



The meaning of this notation will be apparent soon. 

The complete equation describing the dependence on the spatial coordinates is 



P,x dtd> + 



ape 



-di 



a P, x did) + 



di(Dc 



mi 



o?V^ m 4 
M 4 



(3) 



(4) 



(5) 



Instead of solving these equations, the authors of [121 131 US] analyzed (though not solved) 
equation (j3J) ignoring the higher derivative term d t (<p + 3H<p)/m 2 , assuming that it is small. Then 
they analyzed equation applying it to perturbations, ignoring the term d t (Dd)) / rrig, but keeping 
the term di(\3(f))/m 2 , assuming that it is large. Our goal is to see what happens if one performs 
an investigation in a self-consistent way. 

In order to do this, let us temporarily assume that the higher derivative term is absent, which 
corresponds to the limit m g — > oo. In this case our equation for reduces to the equation used in 

pa m na coo 

d t [a 3 P, x d)} = -a 3 V, (P m 4 /M 4 . (6) 
One of the Einstein equations, in the same approximation, is 



H 



- l -{e + p) = -M 4 P, X X = -M A X{X - 1/2) 



(7) 



where e is the energy density and p is the pressure. (We are using the system of units where 
Ml = (87TG)- 1 = 1.) 

The null energy condition (NEC) requires that e +p > 0, and H < 0. Therefore a collapsing 
universe with H < cannot bounce back unless NEC is violated. It implies that the bounce can 
be possible only if P,x becomes negative, P,x < 0, i.e. the field X should become smaller than 
1/2. 

It is convenient to represent the general solution for d)(t) as 



d)(t) = -m 2 t + 7r (t) + ir(xi, t) , 



where vr (t) satisfies equation 



7r + 3Hn 



m 



M 4 



(8) 



(9) 



In this case one can show that the perturbations of the field ir(xi,t) have the following spectrum 
at small values of P,x- 



to 



Px^ . 



(10) 
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This means that P,x plays in this equation the same role as the square of the speed of sound. For 
small P,x, one has 



c 2 = p, x . nr. 



The ghost condensate point P,x = 0, which separates the region where NEC is satisfied and the 
region where it is violated, is the point where the perturbations are frozen. The real disaster 
happens when one crosses this border and goes to the region with P,x < 0, which corresponds 
to < 0. In this area the NEC is violated, and, simultaneously, perturbations start growing 
exponentially, 

7T k (t) ~ exp( VlcfT \k\t)~ exp(v/|P, x | \k\ t) . (12) 

This is a disastrous gradient instability, which is much worse than the usual tachyonic instability. 
The tachyonic instability develops as exp(\/m 2 — k 2 t), so its rate is limited by the tachyonic mass, 
and it occurs only for k 2 < m 2 . Meanwhile the instability (I12p occurs at all momenta k, and the 
rate of its development exponentially grows with the growth of k. This makes it abundantly clear 
how dangerous it is to violate the null energy condition. 

That is why it was necessary to add higher derivative terms of the type of — | (jw^ to the 

ghost condensate Lagrangian [16]: The hope was that such terms could provide at least some 
partial protection by changing the dispersion relation. 

Since we are interested mostly in the high frequency effects corresponding to the rapidly de- 
veloping instability, let us ignore for a while the gravitational effects, which can be achieved by 
taking a(t) = 1, H = 0. In this case, the effective Lagrangian for perturbations 7r of the field <fi in 
a vicinity of the minimum of P(X) (i.e. for small \P,x I) is 



M 4 

■ , ■ . - ■ ■ > (13) 



m 4 



I* 2 - ¥<* < v *> 2 - 2^f w 



The equation of motion for the field tt is 

TT - P, X V 2 7T + — ^dV = . (14) 

rrig 

At small frequencies u, which is the case analyzed in [16], the dispersion relation corresponding 
to this equation looks as follows: 

k 4 

,2 



u 2 = P, x k 2 + ^ . (15) 



m g 



This equation implies that the instability occurs only in some limited range of momenta k, which 
can be made small if the parameter m g is sufficiently small and, therefore, the higher derivative 
therm is sufficiently large. This is the one of the main assumptions of the new ekpyrotic scenario: 
If the violation of the NEC occurs only during a limited time near the bounce from the singularity, 
one can suppress the instability by adding a sufficiently large term — ^ziOir) 2 . (This term must 

have negative sign, because otherwise it does not protect us from the gradient instability. This 
will be important for the discussion in Appendix.) 

Note that one cannot simply add the higher derivative term and take it into account only up 
to some cut-off u 2 ,k 2 < A 2 . For example, if we "turn on" this term only at k 2 < A 2 , it is not 
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going to save us from the gradient instability which occurs at uj 2 = P,x k 2 for all unlimitedly large 
k in the region where the NEC is violated and P,x < 0. 

There are several different problems associated with this scenario. First of all, in order to 
tame the instability during the bounce one should add a sufficiently large term —^(jp) , which 
leads to the emergence of the term ^ziOn) 2 in the equation for 7r. But if this term is large, then 
one should not discard it in the equations for the homogeneous scalar field and in the Einstein 
equations, as it was done in [12], [131 El H5]- 

The second problem is associated with the way the higher derivative terms were treated in 
[TBI [T2| [T3"t dU [T5]. The dispersion relation studied there was incomplete. The full dispersion 
relation for the perturbations in the theory fflBl . (fT4j) is 

uj =P, x k H - . (16) 

m A g 

This equation coincides with eq. (TT5l) in the limit of small u studied in [TBI [121 El El E] • However, 
this equation has two different branches of solutions, which we will present, for simplicity, for the 
case P,x = corresponding to the minimum of the ghost condensate potential P{X): 

to = ± uJi , i — 1, 2, (17) 

where 

wi = \ {^Jml + Ak 2 - m g ) , (18) 
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to 2 = -^m 2 g + Ak 2 + m g j . (19) 

At high momenta, for k 2 ^> m 2 , the spectrum for all 4 solutions is nearly the same 

u w ±|fc| . (20) 

At small momenta, for k 2 <ti m 2 , one has two types of solutions: The lower frequency solution, 
which was found in [16], is 

oj = ± k 2 /m g . (21) 
But there is also another, higher frequency solution, 

u = ± m g . (22) 

The reason for the existence of an additional branch of solutions is very simple. Equation for 
the field <fi in the presence of the term with the higher derivatives is of the fourth order. To specify 
its solutions it is not sufficient to know the initial conditions for the field and its first derivative, 
one must know also the initial conditions for the second and the third derivatives. As a result, a 
single equation describes two different degrees of freedom. 

To find a proper interpretation of these degrees of freedom, one must perform their quanti- 
zation. This will be done in the next two sections. As we will show in these sections, the lower 
frequency solution corresponds to normal particles with positive energy uj = +Ui(k), whereas 
the higher frequency solution corresponds to ekpyrotic ghosts with negative energy —u>2{k). The 
quantity — m g has the meaning of the ghost mass: it is given by the energy uj = —u 2 (k) at k = 
and it is negative. 
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3 Hamiltonian quantization 



We see that our equations for lu have two sets of solutions, corresponding to states with positive 
and negative energy. As we will see now, some of them correspond to normal particles, some of 
them are ghosts. We will find below that the Hamiltonian based on the classical Lagrangian in 
eq. (USD is 

—-3 \u) X a\a k - u 2 c[c k J . (23) 

The expressions for u)\ and io 2 will be presented below for the case of generic c 2 , for c 2 = 
they are given in eqs. (fTBj) and f|T9|) . Both uj\ and uo 2 are positive, therefore a\ and a k are 
creation/annihilation operators of normal particles whereas c k and c k are creation/annihilation 
operators of ghosts. 

We will perform the quantization starting with the Lagrangian in eq. fflBl . with an arbitrary 
speed of sound, c 2 = P,x- The case c s = 1 is the Lorentz invariant Lagrangian. The case 
c 2 s = P,x = is the case considered in the previous section and appropriate to the ghost condensate 
and the new ekpyrotic scenario at the minimum of P{X). 

By rescaling the field tt — > ^tt we have 

L = \ (d^it + (c 2 s - 1) ttAtt - -\(On) 2 ) . (24) 
2 \ rrig J 

This is the no-gravity theory considered in the previous section. Note that the ghost condensate 
set up is already build in, the negative kinetic term for the original ghost is eliminated by the 
condensate. The existence of higher derivatives was only considered in [12], [13|, [HJ [15l [16] as 
a 'cure' for the problem of stabilizing the system after the original ghost condensation. As we 
argued in the previous section, this 'cure' brings in a new ghost, which remained unnoticed in 
[T2"l [TBI dH [TSJ, [12]. In this section, as well as in the next one, we will present a detailed derivation 
of this result. Because of the presence of higher derivatives in the Lagrangian, the Hamiltonian 
quantization of this theory is somewhat nontrivial. It can be performed by the method invented 
by Ostrogradski [23] . 

Thus we start with the rescaled eq. ([TBI 

L = i 
2 

The equation of motion for the field 71 is 



7r 2 -^(V,vr) 2 -^(D7r) 
Trig 



(25) 



n - c 2 V 2 vr + -^n 2 7T = . (26) 

If the Lagrangian depends on the field it and on its first and second time derivatives, the general 
procedure is the following. Starting with L = L(n, 7r, tt), one defines 2 canonical degrees of freedom, 
(qi,Pi) and (q 2 ,P2)- 

dL d dL 
qi = tt , px = — - — — , 
oqi dt oqi 

P2 = — . (27) 
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The canonical Hamiltonian is 

H = piqi + p 2 q2 ~ L(qi, <ti, qi, <?2 
The canonical Hamiltonian equations of motion, 

dH . dH 



dpi 



Pi 



dqi 



i = l,2 



(28) 



(29) 



are standard; they exactly reproduce the Lagrangian equation of motion (1261) . The quantization 
procedure requires promoting the Poisson brackets to commutators which allows to identify the 
spectrum. There are many known examples of the Ostrogradski procedure of derivation of the 
canonical Hamiltonian, see for example [2H [25] . 

The Hamiltonian density constructed by the Ostrogradski procedure for the Lagrangian (J25]) 

is 

H d (x,t) = \\p\ - (p! - q 2 f - m] (p 2 - -^V 2 xgi ) + c 2 s {V x q x f + \ (V* gi) 2 ] . (30) 

The next step in quantization is to consider the ansatz for the solution of classical equations of 
motion in the form 



qt(x, t) 



d 3 k 
(2tt) s 



fl f2 

J k —ik\x I J k a ik2X 



+ 



+ CC 



(31) 



where h\ x = oj\{k) t — k x, and k 2 x = 0J 2 {k) t — k x. We impose the Poisson brackets 

{q i} p j } = 6 ij (32) 

and promote them to commutators of the type 

[q i (x,t),p j (x',t)} = i5 ij 5 3 (x- X 1 ) . (33) 

This quantization condition requires to promote the solution of the classical equation (1311) to the 
quantum operator form, where 



fl = o.k 



P 2 _ _ m 9 



fk = c k 



UJn — CUT 



(34) 



LUn — CUT 



and we impose normal commutation relation both on particles with creation and annihilation 
operators a' and a and ghosts, and c: 



[a k , 4,] = (2ixf5\k - k') , [c k , c\] = (2n) 3 5 3 (k - k') 



Here 



u)i(k ;m g ,c s ) = \k + 



2 , m 9 



k 2 m 2 {1 _ C 2) + ^ 



1/2 



(35) 



(36) 
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and 



2 / T\ V2 

771 I 771 \ 

u 2 (k 2 ; m g , c 2 s )=\k 2 + ^ + \jk*m? g {l - <*) + -f\ (37) 



The Hamiltonian operator acquires a very simple form 

1 f d 3 k / \ f d 3 k / \ 

H quant = -J (wi(4 a * + °*4) ~ W 2 (4 C * + C *4)J = J 7^3 (^4^ ~ ^4 C & + C j 



(38) 

Here the infinite term C is equal to — co>2)(5 3 (0). It represents the infinite shift of the 

vacuum energy due to the sum of all modes of the zero-point energies and is usually neglected 
in quantum field theory. Apart from this infinite c-number this is an expression promised in eq. 

(EH}. 

We now define the vacuum state |0) as the state which is annihilated both by the particle as 
well as by the ghosts annihilation operators, Ofc|0) = c^|0) = 0. Thus the energy operator acting 
on a state of a particle has positive value and on a state of a ghost has the negative valu^j] 

H quant 4|0) = toi{k) 4|0) , H quant 4|0) = -u 2 (k) cj|0> . (39) 

This confirms the physical picture outlined at the end of the previous section. 



4 Lagrangian quantization 

The advantage of the Hamiltonian method is that it gives an unambiguous definition of the 
quantum-mechanical energy operator, which is negative for ghosts. This is most important for our 
subsequent analysis of the vacuum instability in the new ekpyrotic scenario. However, it is also 
quite instructive to explain the existence of the ghost field in new ekpyrotic scenario using the 
Lagrangian approach. The Lagrangian formulation is very convenient for coupling of the model 
to gravity. 

Using Lagrangian multiplier, one can rewrite eq. (|24|) as 



1 / B 2 \ 
L = -[ daird** + (c 2 s - 1) ttAtt + A (B - □ 

2 \ m z J 



tt) . (40) 

Variation with respect to B gives A = — y. After substituting A in fT40l) and skipping total derivative 
we obtain 

L = \ (d^n + (c 2 s - 1) tt Att + ^]+ \d^B d»n 
2 \ ' nig J m l g 



1 Onc could use an alternative way of ghost quantization, by changing the sign of their commutator relations. In 
this case the ghosts would have positive energy, but this would occur at the expense of introducing a nonsensical 
notion of negative probabilities. 
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Introducing new variables a, £ according to 



2B 

CT + 4 = 7T H 



m 3 



cr - £ = 7T , 

and substituting 7r = a — £, £> = D7T = m 2 £ in (jUJ) we obtain 



L k=0 = ^(a 2 -i 2 + m 2 g e) ■ (44) 



L = i (S^o- - + m 2 £ 2 ) + i (c 2 - 1) (a - £)A(a - £) . (42) 

In the case c 2 = 1 we have two decoupled scalar fields: massive with negative kinetic energy and 
massless with positive kinetic energy. 

L cl=1 = l - (d.ad^cr - d„£&>t + m 2 g e) • (43) 
For the homogeneous field (k = mode) the Lagrangian does not depend on c 2 and is reduced 

to 

1 

2 

The relation between the k = mode of the original field n, the normal field a and the ghost field 

a = 7r + £, e=-^- (45) 

When c 2 7^ 1 and 7^ these fields still couple. To diagonalize the Lagrangian in eq. fl24"j) and 
decouple the oscillators we have to go to normal coordinates, similar to the case of the classical 
mechanics of coupled harmonic oscillators. For that we need to solve the eigenvalue problem and 
find the eigenfrequencies of the oscillators. Let us consider the modes with the wavenumbers k. 
For such modes we can perform the following change of variables 

7f k + ^|7r k vr k + u; 2 7r k 

= , , 4k = , , (46) 

m g V ^2 - m 9 a/o;| - wf 

where u>i,a>2 are defined in eqs. fl36l) . fl37j) . In the special case of c 2 = the answers for ui,u)2 
simplify and are shown in eqs. ( fl~8l) and ( fl9|) . After a change of variables we find for these modes 
in the momentum space 

Lcl = 2 (o"k^_k - ^1 o"kO"-k - £k£-k + ^2 fkf-k) • (47) 

The modes of a are normal, and the modes of £ are ghosts. Using this Lagrangian, one can 
easily confirm the final result of the Hamiltonian quantization given in the previous section^ 



2 After we finished this paper, we learned that the Lagrangian quantization of the ghost condensate scenario 
was earlier performed by Aref'eva and Volovich for the case c 2 s = P,x = [26], and they also concluded that this 
scenario suffers from the existence of ghosts. When our works overlap, our results agree with each other. We use 
the Lagrangian approach mainly to have an alternative derivation of the results of the Hamiltonian quantization. 
The Hamiltonian approach clearly establishes the energy operator and the sign of its eigenvalues, which is necessary 
to have an unambiguous proof that the energy of the ghosts is indeed negative and the ghosts do not disappear at 
non-vanishing c 2 s = P,x 7^ when the ekpyrotic universe is out of the ghost condensate minimum. 
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The classical mode <Jk is associated with creation/annihilation operators a k , au of normal particles 
after quantization and the classical mode £k is associated with creation/annihilation operators of 
ghosts particles c\,Ck after quantization. A quantization of the theory in eq. (|4"T|) leads to the 
Hamiltonian in eq. (123]) . 



5 Energy- momentum tensor and equations of motion 

First, we will compute the energy-momentum tensor (EMT) of the Lagrangian (I42p using the 
Noether procedure: 



where 77^ is Minkowski metric. We find 



d^ad u a - d^dvt + ^(c 2 . - l)d\a - £)d u {a - £) 



(48 



The energy density is 



e = Too = -[a 2 + (d lC r) 2 - £ 2 ~ (^) 2 - ^ + (c 2 s ~ 1) - £)) 2 ]. 

For the homogeneous field (k = mode), the energy density can be split into two parts, i.e. a 
normal field part and an ekpyrotic ghost field part: 

e = e a + ££, 



where 



-a 2 > 0, 



1 • 1 



e e = — P - -mi? < . 

Thus the energy of the ghost field £ is negative. 

Up to now we have turned off gravity. In the presence of gravity, the energy-momentum tensor 
of the full Lagrangian (1) in Sec. 2 is calculated by varying the action with respect to the metric: 

. M 'P(X) - <5f + m - d " {md "'- 



T ^ U 9 » v ' 2M' 2 ' ' VYV M' 2 

+M 4 m~ 4 P, x d^d v (j) + M'- 2 (^(D0) d v (j> + d v {U(j>)d^ 
M' 2 Y 2 



9fj,u 



-M 4 P(X) 



+ V(<f>) -d Q Yd° 



+M 4 m~ 4 P,x d^d u <P + d^Y d v (f> + d v Y d : 



where 



y = M'- 2 n 



(49) 
(50) 
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From this, for a homogeneous, spatially flat FRW space time we have the energy density 

M' 2 Y 2 ■ ■ 

e = M 4 (2P, x X-P(X)) + V((f)) — + YCI) (51) 



M' 2 Y 2 



and the pressure 

p = M*P(X)-V(o)l —7—- Yo. (52) 

so that 



H = ~(e + p) = -M 4 P, X X-Y'<p. (53) 

Note that in the homogeneous case in the absence of gravity the ekpyrotic ghost field £ as 
defined in eq. (|45|) is directly proportional to the field Y: 

C = ^ Y . (54) 
M 2 K J 

The closed equations of motion, which we used for our numerical analysis, are obtained as 
follows: 

<t>(P, x +2XP, xx ) + 3HP lX <j ) +—(Y + 3HY) = -V^m^/M 4 , 

H = -M 4 P, X ^-Y'<j>, (55) 
M'- 2 (4> + 3H<j)) — Y . 

Here X = <p 2 /2m 4 . 

In these equations the higher derivative corrections appear in the terms containing the deriva- 
tives of Y. The last of these equations shows that Y — >• in the limit M' — > oo (i.e. m g — > oo), 
and then the dynamics reduces to one with no higher derivative corrections. 

The closed equations of motion for it coupled to gravity are obtained by expanding (155]) and 
linearizing with respect to ir and Y. Then we get 

4 

n + 3Hir + — (Y + 3HY) = -V,^m 4 /M 4 , 

7f + 3#7T = M' 2 Y + 3Hm 2 , (56) 

H = —— + Ym 2 . 

2m 2 
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6 On reality of the bounce and reality of ghosts 



Using the equations derived above, we performed an analytical and numerical investigation of the 
possibility of the bounce in the new ekpyrotic scenario. We will not present all of the details of this 
investigation here since it contains a lot of material which may distract the reader from the main 
conclusion of our paper, discussed in the next section: Because of the existence of the ghosts, this 
theory suffers from a catastrophic vacuum instability If this is correct, any analysis of classical 
dynamics has very limited significance. However, we will briefly discuss our main findings here, 
just to compare them with the expectations expressed in [121 131 HH US] • 

Our investigation was based on the particular scenario discussed in [13], [15] because no explicit 
form of the full ekpyrotic potential was presented in [121 EH] • The authors of [131 [15] presented 
the full ekpyrotic potential, but they did not fully verify the validity of their scenario, even in the 
absence of the higher derivative terms. 

Before discussing our results taking into account higher derivatives, let us remember several 
constraints on the model parameters which were derived in [TJl [321 EH EES]- We will represent 
these constraints in terms of the ghost condensate mass instead of the parameter M', for K = 1. 
In this case the stability condition (7.19) in [13] (see also [l2l[Tl]) reads: 

!4<-<l«|. (57) 
ti | trig 

It was assumed in [13] that the bounce should occur very quickly, during the time At < |i/ | _1 ~ 



l/y/p\Vmin\- Here _H is the Hubble constant at the end of the ekpyrotic state, just before it 
start decreasing during the bounce, p ~ 1CT 2 , and V m i n is the value of the ekpyrotic potential 
in its minimum. During the bounce one can estimate AH ~ \H \ ~ HAt < -M^ because we 

assume, following [13], that At < |-H | _1 ) an d we assume an approximately linear change of H 
from — \Hq\ to |, 
quite restrictive, 



from — \Hq\ to | -Ho I - This means that > \Hq\. In this case the previous inequalities become 



M 4 

\H,\ < — < H ■ (58) 

This set of inequalities requires that the stable bounce is not generic; it can occur only for a 
fine-tuned value of the ghost mass, 

M 4 M 4 

m g ~ ~ — / . (59) 



\fWrr, 



The method of derivation of these conditions required an additional condition to be satisfied, 
I -Ho I ~ \fp\y m in\ M 2 , see Eqs. (8.8) and (8.17) of Ref. [13]. This condition is satisfied for 

m g > M 2 . (60) 



Whereas the condition (159"]) seems necessary in order to avoid the development of the gravita- 
tional instability and the gradient instability during the bounce for m g 3> M 2 , it is not sufficient, 
simply because the very existence of the bounce may require m g to be very much different from 
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Figure 1: The "new ekpyrotic potential," see Fig. 3 in [13] and Fig. 6 in [15]. The cosmological evolution 
in this model results in a universe with a permanently growing rate of expansion after the bounce, which 
is unacceptable. 



its fine-tuned value m n 



Ai 4 



Indeed, our investigation of the cosmological evolution in 



this model shows that generically the bounce does not appear at all, or one encounters a singular 
behavior of (j) because of the vanishing of the term P,x +2XP,xx m dS5D, or one finds an unstable 
bounce, or the bounce ends up with an unlimited growth of the Hubble constant, like in the Big 
Rip scenario [27]. Finding a proper potential leading to a desirable cosmological evolution requires 
a lot of fine-tuning, in addition to the fine-tuning already described in [T3J, [15] . 

For example, the bounce in the model with the "new ekpyrotic potential" described in [T3], [15] 
and shown in Fig. [1] results in a universe with a permanently growing rate of expansion after the 
bounce, which would be absolutely different from our universe. To avoid this disaster, one must 
bend the potential, to make it approaching the value corresponding to the present value of the 
cosmological constant, see Fig. [2j This bending should not be too sharp, and it should not begin 
too early, since otherwise the universe bounces back and ends up in the singularity. Fig. [3] shows 
the bouncing solution in the theory with this potential. 




Figure 2: An improved potential which leads to a bounce followed by a normal cosmological evolution. 
We do not know whether this extremely fine-tuned potential can be derived from any realistic theory. 



Our calculations clearly demonstrate the reality of the ekpyrotic ghosts, see Fig. H[ which 
shows the behavior of the ghost-related field Y = ^£ near the bounce. The oscillations shown in 
Fig. H] represent the ghost matter with negative energy, which was generated during the ekpyrotic 
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Figure 3: The behavior of the Hubble constant H(t) near the bounce, which occurs near t = 18. To verify 
the stability of the universe during the bounce, one would need to perform an additional investigation 
taking into account the ghost field oscillations shown in Fig. HJ 

Y 



4 - 



2 - 



-2 - 



-4 . 




Figure 4: Ekpyrotic ghost field oscillations. 

collapse. We started with initial conditions Y = Y = 0, i.e. in the vacuum without ghosts, and 
yet the ghost-related field Y emerged dynamically. It oscillates with the frequency which is much 
higher than the rate of the change of the average value of the field (p. 

This shows that the ekpyrotic ghost is not just a mathematical construct or a figment of 
imagination, but a real field. We have found that the amplitude of the oscillations of the ghost 
field is very sensitive to the choice of initial conditions; it may be negligibly small or very large. 
Therefore in the investigation of the cosmological dynamics one should not simply consider the 
universe filled with scalar fields or scalar particles. The universe generically will contain normal 
particles and ghost particles and fields with negative energy. The ghost particles will interact with 
normal particles in a very unusual way: particles and ghosts will run after each other with ever 
growing speed. This regime is possible because when the normal particles gain energy, the ghosts 
loose energy, so the acceleration regime is consistent with energy conservation. This unusual 
instability, which is very similar to the process to be considered in the next section, can make it 
especially difficult to solve the homogeneity problem in this scenario. 
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7 Ghosts, singularity and vacuum instability 



It was not the goal of the previous section to prove that the ghosts do not allow one to solve the 
singularity problem. They may or may not spoil the bounce in the new ekpyrotic scenario. How- 
ever, in general, if one is allowed to introduce ghosts, then the solution of the singularity problem 
becomes nearly trivial, and it does not require the ekpyrotic scenario or the ghost condensate. 

Indeed, let us consider a simple model describing a flat collapsing universe which contains a 
dust of heavy non-relativistic particles with initial energy density pm, and a gas of ultra-relativistic 
ghosts with initial energy density — p g < 0. Suppose that at the initial moment t = 0, when the 
scale factor of the universe was equal to a(0) = 1, the energy density was dominated by energy 
density of normal particles, pu — Pg > 0. The absolute value of the ghost energy density in the 
collapsing universe grows faster than the energy of the non-relativistic matter. The Friedmann 
equation describing a collapsing universe is 



In the beginning of the cosmological evolution, the universe is collapsing, but when the scale 
factor shrinks to cibounce = ^j, the Hubble constant vanishes, and the universe bounces back, thus 
avoiding the singularity. 

Thus nothing can be easier than solving the singularity problem once we invoke ghosts to help 
us in this endeavor, unless we are worried about the gravitational instability problem mentioned 
in the previous section. Other examples of the situations when ghosts save us from the singularity 
can be found, e.g. in [28], where the authors not only study a way to avoid the singularity with 
the help of ghosts, but even investigate the evolution of metric perturbations during the bounce. 
So what can be wrong with it? 

A long time ago, an obvious answer would be that theories with ghosts lead to negative 
probabilities, violate unitarity and therefore do not make any sense whatsoever. Later on, it was 
realized that if one treats ghosts as particles with negative energy, then problems with unitarity 
are replaced by the problem of vacuum stability due to interactions between ghosts and normal 
particles with positive energy, see, e.g. [23 E01 ED E21 E31 EH E5]- Indeed, unless the ghosts are 
hidden in another universe [SU] , nothing can forbid creation of pairs of ghosts and normal particles 
under the condition that their total momentum and energy vanish. Since the total energy of ghosts 
is negative, this condition is easy to satisfy. 

There are many channels of vacuum decay; the simplest and absolutely unavoidable one is due 
to the universal gravitational interaction between ghosts and all other particles, e.g. photons. An 
example of this interaction was considered in [32J, see Fig. |5j Nothing can forbid this process 
because it does not require any energy input: the positive energy of normal particles can be 
compensated by the negative energy of ghosts. 

An investigation of the rate of the vacuum decay in this process leads to a double-divergent 
result. First of all, there is a power-law divergence because nothing forbids creation of particles 
with indefinitely large energy. In addition, there is also a quadratic divergence in the integral over 
velocity [321 El]- This leads to a catastrophic vacuum decay. 




(61) 
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Figure 5: Vacuum decay with production of ghosts £ and usual particles 7 interacting with each 
other by the graviton exchange. 

Of course, one can always argue that such processes are impossible or suppressed because of 
some kind of cutoff in momentum space, or further corrections, or non-local interactions. However, 
the necessity of introducing such a cut-off, or additional corrections to corrections, after introducing 
the higher derivative terms which were supposed to work as a cutoff in the first place, adds a lot 
to the already very high price of proposing an alternative to inflation: First it was the ekpyrotic 
theory, then the ghost condensate and curvatons, and finally - ekpyrotic ghosts with negative 
energy which lead to a catastrophic vacuum instability. And if we are ready to introduce an 
ultraviolet cutoff in momentum space, which corresponds to a small-scale cutoff in space-time, 
then why would we even worry about the singularity problem, which is supposed to occur on an 
infinitesimally small space-time scale? 

In fact, this problem was already emphasized by the authors of the new ekpyrotic scenario, 
who wrote [13J: 

"But ghosts have disastrous consequences for the viability of the theory. In order to regulate 
the rate of vacuum decay one must invoke explicit Lorentz breaking at some low scale [32]. In any 
case there is no sense in which a theory with ghosts can be thought as an effective theory, since 
the ghost instability is present all the way to the UV cut-off of the theory." 

We have nothing to add to this characterization of their own model. 

Acknowledgments: It is a pleasure to thank B. Craps, P. Creminelli, G. Horowitz, N. Kaloper, 
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discussions. The work by R.K. and A.L. was supported in part by NSF grant PHY-0244728 and 
by the Alexander-von-Humboldt Foundation. The work by J. U.K. was supported by the German 
Academic Exchange Service. 
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8 Appendix. Exorcising ghosts? 



After this paper was submitted, one of the authors of the new ekpyrotic scenario argued [36J that, 
according to [37J , ghosts can be removed by field redefinitions and adding other degrees of freedom 
in the effective UV theory [37]. Let us reproduce this argument and explain why it does not apply 
to the ghost condensate theory and to the new ekpyrotic scenario. 

Refs. [361 [37] considered a normal massless scalar field with Lagrangian density in (—,+,+, +) 
signature o 

£ = -I(90) 2 + ^W) 2 - V int (<P) , (62) 

where a = ±1, and Vi nt is a self-interaction term. This theory is similar to the ghost conden- 
sate/new ekpyrotic theory in the case a = — 1, c s = 1, see eqs. (1) and (12^1) . The sign of a 
is crucially important: the term (D0) 2 would not protect this theory against the gradient 

instability in the region with the NEC violation. 

Note that in notation of [361 E2]; m g = A, which could suggest that the ghost mass is a 
UV cut-off, and therefore there are no dangerous excitations with energies and momenta higher 
than m g . However, this interpretation of the theory fl62l) would be misleading. Upon a correct 
quantization, this theory can be represented as a theory of two fields without the higher derivative 
non-renormalizable term ^2 (He/)) 2 , see Eq. f l42l) . One can introduce the UV cut-off A when 
regularizing Feynman diagrams in this theory, but there is absolutely no reason to identify it 
with m g \ in fact, the UV cut-off which appears in the regularization procedure is supposed to be 
arbitrarily large, so the perturbations with frequencies greater than m g should not be forbidden. 

Moreover, as we already explained in Section 2, one cannot take the higher derivative term 
into account only up to some cut-off uj 2 ,k 2 < A 2 . If, for example, we "turn on" this term only at 
k 2 < A 2 , it is not going to protect us from the gradient instability, which occurs at u 2 = P,x k 2 
for all indefinitely large k in the region where the NEC is violated and P,x < 0. Note that this 
instability grows stronger for greater values of momenta k. Therefore if one wants to prove that the 
new ekpyrotic scenario does not lead to instabilities, one must verify it for all values of momenta. 
Checking it for uj 2 ,k 2 < m 2 is insufficient. Our results imply that if one investigates this model 
exactly in the way it is written now (i.e. with the term — 7^2- (D0) 2 ), it does suffer from vacuum 
instability, and if we discard the higher derivative term at momenta greater that some cut-off, the 
instability becomes even worse. Is there any other way to save the new ekpyrotic scenario? 

One could argue [361 137] that the term (D^) 2 is just the first term in a sum of many higher 
derivative terms in an effective theory, which can be obtained by integration of high energy degrees 
of freedom of some extended physically consistent theory. In other words, one may conjecture that 
the theory can be made UV complete, and after that the problem with ghosts disappears. However, 
not every theory with higher derivatives can be UV completed. In particular, the possibility to 
do it may depend on the sign of the higher derivative term [19] . 

According to [37] , the theory (1621) is plagued by ghosts independently of the sign of the higher 

3 In our paper we used the signature (+, — , — , — ), so some care should be taken when comparing the equations. 
Note that this does not change the sign of the higher derivative term ~ a(D0) 2 ; the ghost condensate/ekpyrotic 
theory corresponds to a = — 1. 
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derivative term in the Lagrangian. One can show it by introducing an auxiliary scalar field x an d 
a new Lagrangian 

C' = -\(d<P) 2 -ad^<p-\am 2 gX 2 -V,M), (63) 
which reduces exactly to C once x is integrated out. C is diagonalized by the substitution 
= 0'- a-X- 

C! = -|(<90') 2 + \{d X f - \a m 2 g x 2 ~ VUt, X) , (64) 

which clearly signals the presence of a ghost: x has a wrong-sign kinetic term. 

Then the authors of [37] identified x as a tachyon for a = —1, suggesting that in this case x 
has exponentially growing modes. However, this is not the case: due to the opposite sign of the 
kinetic term for the Y-field, the tachyon is at a = +1, not at a = — 1. Indeed, because of the flip 
of the sign of the kinetic term for the field x-, hs equation of motion has a solution x ~ e ±l ^ t - kx ) 
with 

- (to 2 -k 2 )=am 2 g . (65) 

For the field with the normal sign of the kinetic term, the negative mass squared would mean 
exponentially growing modes. But the flip of the sign of the kinetic term performed together 
with the flip of the sign of the mass term does not lead to exponentially growing modes [25J ED]- 
Based on the misidentification of the negative mass of the field with the wrong kinetic terms as 
a tachyon, the authors choose to continue with the a = +1 case in eq. fl62l) . Starting from this 
point, their arguments are no longer related to the ghost condensate theory and the new ekpyrotic 
theory, where a = — 1. We will return to the case a = — 1 shortly. 

For the a = 1 case they argued that the situation is not as bad as it could seem. They proposed 
to use the scalar field theory eq. (I6"2"|) at energies below m g , and postulated that some new degree of 
freedom enters at k > m g and takes care of the ghost instability. The authors describe this effect 
by adding a term — (dx) 2 to construct the high energy Lagrangian. For V\ nt = they postulate 

C$ =£- (d X ) 2 = -i(<90) 2 - dtfPfj, - {dx? -\m\ x 2 (66) 

and use the shift = — x to get a simple form of a UV theory. This trick reverses the sign of 
the kinetic term of the field x-, an d the ghost magically converts into a perfectly healthy scalar 
with mass m g : 

C a ^ = -\(d^) 2 -\(dx) 2 -\m 2 g x 2 ■ (67) 

One may question validity of this procedure, but let us try to justify it by looking at the 
final result. Consider equations of motion for x from eq. (1661) and solve them by iteration in the 
approximation when □ in 2 - 
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□0 Dx 06 2 6 , . 

— "-) _ iiyj _ ivy) _ ivy) ri 



m g m g m g m g 



Now replace x m ec i- (|66l) by its expression in terms of OS as given in eq. (|68|) . The result is 
our original Lagrangian (|62|) . plus some additional higher derivative terms, which are small at 
|D| m 2 , i.e. at \u 2 — k 2 \ m 2 . Thus one may conclude that, for a = 1, the theory (|62j) . which 
has tachyonic ghosts, may be interpreted as a low energy approximation of the UV consistent 
theory fl6T|) . 
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Now let us return to the ghost condensate/new ekpyrotic case. To avoid gradient instabilities 
in the ekpyrotic scenario, the sign of the higher derivative term in eq. (|62l) has to be negative, 
a = — 1, see eq. (1) and also eq. (fT5j) and the discussion below it. This means that one should 
start with eq. (J62"j) with a = — 1. 

This theory is not tachyonic, but, as we demonstrated by performing its Hamiltonian quanti- 
zation, it has ghosts, particles with negative energy, in its spectrum. Can we improve the situation 
by the method used above? Let us start with the same formula, CJ — {dx) 2 , as in a = +1 case: 

C a - 1 =£,'- (d X ) 2 = -\(d<pf + d» X dy - (d X ) 2 + \m 2 gX 2 (69) 

and replace x by the iterative solution of its equation of motion 

□0 Dy 06 2 6 
m g m g m g m g 

Thus, up to the terms which are small at |d| <C m 2 , the theory with the Lagrangian ( 1691) does 
reproduce the model (1621) with a = —1, up to higher order corrections in |D|/m^. The theory (1691) 
can be also written as 

^w 1 = -\i^f ~ \idx? + \ m 2 g x 2 ■ (71) 
where 6 = 6 + X- The sign of the kinetic term of both fields is normal, but the mass term still 
has the wrong sign, which leads to the tachyonic instability 5x ~ exp \J m 2 g — k 2 1. Therefore the 
cure for the ghost instability proposed in [36j [37] does not work for the case a = — 1 of the ghost 
condensate/ekpyrotic scenario. 

Moreover, the procedure described above is valid only for \uj 2 — k 2 \ <ti m g . Meanwhile the 
gradient instability of the ekpyrotic theory in the regime of null energy condition violation (c 2 < 0) 
is most dangerous in the limit k 2 —>■ oo, where this procedure does not apply, see ffTST) . This agrees 
with the general negative conclusion of Refs. [191 12Q1 [ZD, (22] with respect to the theories of this 
type. 

In this Appendix we analyzed the Lorentz-invariant theory (162!) because the argument given 
in [3TJI I3T] was formulated in this context. The generalization of our results for the ghost conden- 
sate/new ekpyrotic case is straightforward. Indeed, our results directly follow from the correlation 
between the sign of the higher derivative term in (1631) and the sign of the mass squared term in 
( IMj) . One can easily verify that this correlation is valid independently of the value of c 2 , i.e. at 
all stages of the ghost condensate/new ekpyrotic scenario. 

To conclude, our statement that the ghost condensate theory and the new ekpyrotic scenario 
imply the existence of ghosts is valid for the currently available versions of these theories, as they 
are presented in the literature. In this Appendix we explained why the recent attempts to make the 
theories with higher derivatives physically consistent [361 EZ] do not apply to the ghost condensate 
theory and the new ekpyrotic scenario. One can always hope that one can save the new ekpyrotic 
scenario and provide a UV completion of this theory in some other way, but similarly one can 
always hope that the problem of the cosmological singularity will be solved in some other way. 
Until it is done, one should not claim that the problem is already solved. 
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